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Abstract 

I review the status of CP violation in the Standard Model from the combination of flavour constraints within the 
CKMfitter frequentist approach and I describe studies of New Physics restricted to the AF = 2 sector to explain recent 
results on neutral-meson mixing. All results have been obtained using data available for the Winter 2012 conferences. 



1. Standard Model global fit 

In the Standard Model (SM), the weak charged- 
current transitions mix different quark generations, 
which is encoded in the unitary Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. In the case of three genera- 
tions of quarks, the physical content of this matrix re- 
duces to four real parameters, among which one phase, 
the only source of CP violation in the quark sector. 
These four real parameters are defined in an phase- 
convention independent way: 
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As these parameters have been determined more and 
more precisely over the last decade, the attention has 
shifted from the pure metrology of the Standard Model 
to the investigation of deviations in flavour physics, tell- 
tale signs of New Physics (NP). The CKMfitter group 
follows this programme within a frequentist approach, 
using the Rfit model to describe systematic uncertain- 
ties 0"). A table of the inputs summarising the SM 
global fit is presented in Table [T] Low-energy strong 
interactions constitute a central issue in flavour physics, 
which explains the need for accurate inputs for hadronic 
quantities such as decay constants, form factors, bag pa- 
rameters. . . We mostly rely on Lattice QCD simulations, 
with a specific averaging procedure (Educated Rfit) to 
combine the results from different collaborations. A 
similar approach is followed in order to combine the in- 
clusive and exclusive determinations of \V u b\ and \V c b\- 
An input which has varied recently is the angle y, 



whose measurement is based on the interference be- 
tween the colour-allowed B — > D°K~ and colour- 
suppressed B — > D°K~ decays (as well as modes with 
similar flavour structure). The accuracy of the method is 
driven by the size of the ratio r B = |A SU pp le s S ed I / 1^4 allowed I 
typically of order 0.1-0.2. The different methods try to 
improve on this ratio by different choices of D decay 
channels: GLW (D into CP eigenstates), ADS (D (t) into 
doubly Cabibbo-suppressed states), GGSZ (D w into 3- 
body state and Dalitz analysis). Each method allows 
for a simultaneous determination of y and the relevant 
hadronic quantities, in particular the ratio of amplitudes 
rs and the relative phase between the suppressed and 
allowed amplitudes 6b- Recent measurements for the 
ADS method from CDF (2) and Belle in Summer 
201 1 and from LHCb @ in Winter 2012 together with 
a more sophisticated statistical treatment of the nuisance 
parameters (constrained supremum) have led to a signif- 
icant improvement for the input of the global fit corre- 
sponding to the combination of the three methods. 
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The fit results in the (p, 77) plane are shown in Fig. [T[ 
with the four real parameters describing the CKM ma- 
trix at 68% CL: 
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Even though there is a good overall agreement, the SM 
global fit exhibits a discrepancy between its various ob- 
servables @. Indeed, the global fit Xmi„ drops by 2.8cr 
if sin 2/3 or B — > tv is removed from the inputs. Before 
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Quantity 


Experimental information 


Theoretical input 


\V„d\ 


superallowed B decays 


Ref. 


\Vus\ 


WA K(2,K(2, Ji(2, t — * Kv T , T — > 7TV T 


/ + (0) = 0.963 + 0.003 + 0.005 




WA 


B K = 0.733 ± 0.003 ± 0.036 


\v ub \ 


b -> utv, B -> tt^v 


\Vub\ ■ 10 3 = 3.92 ± 0.09 ± 0.45 


Wcb\ 


ft -> c^v, B -> D ( *^v 


|y ci | . 10 3 = 40.89 ± 0.38 ± 0.59 




WA mixing 


B B JB Bi = 1.024 ± 0.013 ± 0.015 


Am, 


WA Bj-Bj mixing 


B Bs = 1.291 ± 0.025 ± 0.035 


P 


WA J/ifrK™ 




a 


WA 7T7r, p7T, pp 


isospin 


y 


WA B — > DWZ W 


GLW/ADS/GGSZ 


BR(B -> tv) 


(1.73 ±0.35)- 10~ 4 


f B Jf Bd = 1.218 ± 0.008 ± 0.033 






f B , = 229 + 2 + 6 MeV 



Table 1 : Summary of the inputs for the SM global fit. WA stands for "World average" from PDG or HFAG. For further details, see Ref. [T] and 
references therein. 



claiming any effect from New Physics, one should con- 
sider possible sources for such a discrepancy within the 
SM. On the experimental side, both sin 2B or B — » tv 
showed a good agreement between Babar and Belle. 
On the theoretical side, the issue is not restricted to the 
value of f Bj since one observes a 2.9cr discrepancy from 
the ratio BR(B — > Tv)/Amj, indicating that one needs 
to modify f B(t without affecting f Bd ^B Bd . Three solu- 
tions for this discrepancy could be: a) a change in the 
measurement of B — > tvQ b) a correlated change in lat- 
tice estimates of f Bd and B Br with a value of B B j much 
smaller than 1 (vacuum saturation approximation), c) 
the presence of NP affecting the interpretation of the ob- 
servables in terms of the CKM parameters. 



2. B s -> mx 



The dileptonic decay of the B s meson is among the 
most appealing laboratories to study scalar couplings in 
addition to the SM ones, and complements the current 
effort in measuring B — > K(*)/i/j. observables. The ex- 
perimental limits set by the Tevatron and LHC experi- 
ments on the dimuonic branching ratios ||9] [TO] QT) are 
already giving significant constraints on scenarios be- 
yond the Standard Model. A main limitation in the cur- 
rent prediction arises from the knowledge of the decay 



'This is in particular suggested by the recent measurement of 
BR{B -> tv) from Belle presented at ICHEP 2012 (§']. We do not 
include this result in these proceedings which correspond to the re- 
sults available in Spring 2012. We invite the reader to browse through 
Ref. Q to learn more on the impact of this new result. 



constant f Bs , as seen from the master formula: 
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where Y is the Inami-Lim function at NLO in 0(a s ) Q. 
m, is the top quark mass defined in the MS scheme, con- 
sidered at a scale /j t = 0(m t ), Gf is the Fermi constant, 
sin 2 ffff the electroweak mixing angle, and a em the elec- 
tromagnetic coupling constant at the Z pole, where we 
use the tree-level formula [ 12 1: 



aemiMl) = - ^l2G F M 2 w sin 2 6 W = 1/127.46 



(6) 



There are higher-order corrections to this formula, both 
in the electroweak and strong sectors, which we ac- 
count for by varying the renormalisation scale fi t be- 
tween m,/2 and 2m,. The global SM fit yields an ac- 
curate prediction for the branching ratio: 

io- 9 



BR[B S -» aiVIsm = (3.64^) ■ 



(7) 



If we compare this prediction with other values in 
the literature |[T3l [141 Q3), we find predictions span- 
ning from 3.2 to 3.6, which may seem worrying con- 
sidering that LHCb and CMS bounds (4.5 ■ 10~ 9 and 
7.7 ■ 10~ 9 at 95% CL respectively) are getting close to 
these values. A first source of discrepancy, discussed 
in Ref. lfT3l . corresponds to the choice of schemes for 
the electroweak mixing angle and the top quark mass. 
Any scheme can be chosen a priori, but one would like 
to pick up the one minimising the higher-order elec- 
troweak corrections. But the set of NLO electroweak 
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Figure 1: On the left: 95% CL constraints on the unitarity triangle from the SM global fit.. On the right: discrepancy between the measured values 
of [sin 28, Br(B — » tv)] (cross) and the preferred values by the SM global fit (coloured region). 



corrections to this process is not wholly known, so that 
we do not know the preferred scheme, and different 
schemes seem to lead to a significant variation in the 
central values. The case of K — > jivv, where the NLO 
electroweak corrections are known, led the authors of 
Ref. Ifl3l to prefer the scheme used here. However this 
leaves open the question of the systematic uncertainty 
due to our lack of knowledge about electroweak NLO 
corrections, which is not accounted for in the current 
predictions of the branching ratio [including Eq. ^7J], 



In any case, the difference between refs. lfT3lfT4 l and 
our result cannot be ascribed to a difference of schemes. 
But there are also differences in the inputs used in these 
references and in Eq. Q. Actually, when looking at the 
prediction from the SM global fit presented in Eq. (j7]), 
one should take into account the fact that the global fit 
constrains not only the value of the CKM matrix ele- 
ments, but also some of the hadronic inputs. Indeed, 
fg is constrained indirectly by Am s and Z?g i as they are 
both known with a better accuracy from our average of 
lattice QCD results: the respective relative uncertainties 
for our inputs are 6 %, 3 % and 0.2 % for /J , Bb, and 
Am,. This leads the global fit to predict fg s = 234+jj 
MeV more accurately that the input in Table [T] The 
difference between the results of the SM global fit and 
Ref. lPT3ll for BR(B S — > fi/u) can be illustrated by taking 
the ratio p - BR[B S — > /i + /i~] SM /Am s to eliminate f B) 



in favour of Bg and Am s 
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and insert our best-fit values for the inputs, to be com- 
pared with Ref. Ift3l : 



Value 
B 



m t (m,) (GeV) 
Am, (ps _1 ) 

tb s (ps) 
Br(B s -» fifi) 



Winter 2012 CKMfitter 
1.252 
165.1 
17.73 
1.472 
3.6 • 10" 9 



G3 

1.33 
163.5 
17.77 
1.425 
3.2 • 10" 9 



Hence, the choice of inputs in Ref. Ifl3l and our im- 
proved knowledge of hadronic inputs due to the global 
fit yields the span of predictions in the literature. 

Let us add that one should also include the effect 
of AF.5 + to compare the theoretical branching ratio 
(computed at t = 0, i.e., without mixing) and the mea- 
surement at hadronic machines (integrating over time 
and thus including mixing) lfl6l [T71 [T8l . In the case of 
the Standard Model, the former should be 9% below the 
latter, bringing the LHCb upper limit even closer to the 
SM value predicted from the global fit Eq. (j7]). 

3. New Physics in AF = 2 operators: Mu 

In addition to the discrepancy between B — > tv and 
sin2y6 present at the time of the Winter 2012 update, 
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a large set of information has been obtained concerning 
B S B S mixing, as well as intriguing results for the dimuon 
charge asymmetry. It is interesting in this context to 
look for NP in mixing, i.e., AB — 2 operators, assuming 
that tree decays are not affected by NP effects |[T9ll20ll . 
BqB q oscillations (with q — d or q — s) are described 
by a Schrodinger equation with an evolution matrix for 
\B q (t)), \B q (t)} of the form M q -iT q /2, with the hermitian 
mass and decay matrices M q and Y q . 

The physical eigenstates \B q H ) and with masses 
M q H , M q L and decay rates F^, are obtained by diago- 
nalizing M q - iT q /2. The description of the B q B q oscil- 
lations involves the three physical quantities \M q 2 \, \F q 2 \ 
and the CP phase <p q = arg(-Mj 2 /Fj 2 ). The average B q 
mass and width is M Bq and T Bi , and the mass and width 
differences between B q L and B^ H can be written as 

AM q = M% -M 9 L = 2\M\ 2 \, 
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A third quantity probing mixing is the semileptonic CP 
asymmetry, 



2 1 



Im 
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M q 



K 2 I 



sin< 



a q s is the CP asymmetry in flavour-specific B q 
cays (/ being chosen so that the decays B q - 



(11) 

> / de- 
/ and 

B q — > / are forbidden), denoted a 1 ^ for / semileptonic. 

These quantities are expected to be affected by NP 
differently. While M q n coming from box diagrams is 
very sensitive to NP both for Bj and B s , stems from 
tree-level decays and possible NP effects are generally 
expected to be smaller than the hadronic uncertainties 
from decay constants and bag parameters. Assuming 
that NP does not enter tree-level decays (a more spe- 
cific definition can be found in Ref. [19]), a first natural 
model-independent approach to NP in AF = 2 consists 
in changing the magnitude and/or the phase of M q 2 only. 
It is convenient to define the NP complex parameters A q 
and ^ (q = d, s) as M q 2 = M™' 9 • A q , with 



A, = \Ag\e 1 * 



a J.SM , ,A 

b q = <p q +4> q . 



(12) 



New Physics in M q 2 will not only affect the neutral- 
meson mixing parameters, but also the time-dependent 
analyses of decays corresponding to an interference be- 
tween mixing and decay. 

The following inputs of the SM global fit are consid- 
ered to be free from NP contributions in their extraction 
from experiment: \V ud \, \V US \, \V ub \, \V cb \ and y. We also 
assume that the leptonic decay B — > rv is SM-like. Un- 
der these assumptions, a reference unitarity triangle can 



be constructed, with two solutions for the apex of the 
unitarity triangle (corresponding to the usual solution 
and the symmetric one with respect to the origin of the 
(p, rj) plane). On the other hand, we use several observ- 
ables, affected by NP in mixing, to determine Ad, A s : 
the oscillation frequencies Amj, Am s , the lifetime dif- 
ference ATd, the time dependent asymmetries related to 
4>B d ,tpB,, the asymmetries «g L , aL, A$l, and finally a 
(from interference between decay and mixing). 

In summer 2010 the CDF and D0 analyses of B s — > 
J/i//(f> pointed towards a large negative value of <p^ , while 
simultaneously being consistent with the SM due to 
large errors. A large (pf < could accommodate D0 
large negative value for the semileptonic CP asymme- 
try reading A$ L = 0.6«g L + OAaL GTIl in terms of 
the individual semileptonic CP asymmetries in the Bd 
and B s systems. Moreover, the discrepancy between 
B(B — > rv) and the mixing-induced CP asymmetry in 
Bd J/iffK can be removed with <p^ < 0. The al- 
lowed range for <f>^ implies a contribution to A$ L with 
the right (i.e. negative) sign. In our 2010 analysis |[T9l 
we have performed a simultaneous fit to CKM parame- 
ters, A s and in three generic scenarios in which NP is 
restricted to AF - 2 flavour-changing neutral currents. 
The most general scenario (called Scenario I) treated A s 
and A^ independently, corresponding to NP with arbi- 
trary flavour structure. At that time, we found an excel- 
lent fit in this scenario, with all discrepancies alleviated 
through A^ j + 1. The recent LHCb measurement 11221 
of the CP phase ^ from A™ X (B. S J/ip(f>) does not 
permit large deviations of <f>^ from zero anymore, as in- 
dicated also by CDF [23 1. In Ref. [20|, we have recon- 
sidered our NP scenarios with these additional data. 

In Fig. [2j we show the results in the complex A^ 
and A, planes. Ad and A s are taken as independent 
in this scenario, but some of the constraints correlate 
them, such that Asl from the inclusive dimuon asym- 
metry, and the ratio Am^/Am,. The figures should be 
understood as two-dimensional projections of a single 
multidimensional fit, and not as independent computa- 
tions. On the left panel, The constraint from Anid in the 
ReAj - ImAd plane shows two rings, related to the two 
allowed solutions for the reference triangle in the p - fj 
plane. In the combined fit, the inner ring in the complex 
Ad plane is disfavoured, leading to an allowed region 
for \Ad\ compatible with the SM value A^ = 1 at the 3cr 
level. The NP phase 0*, mainly driven by the discrep- 
ancy between BR(B — > rv) and sin 2B, shows a clear de- 
viation from the SM at the moment of the Winter 2012 
update - a fit without BR(B — > tv) (or with an updated 
value following Ref. [8Q yields (p& d in agreement with 
the SM. On the right panel, we see that the recent LHCb 
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Figure 2: Constraint on the complex parameter Aj and A s . For the individual constraints the coloured areas represent regions with CL < 68.3 %. 
For the combined fit the red area shows the region with CL < 68.3 % while the two additional contour line inscribe the regions with CL < 95.45 %, 
and CL < 99.73 %, respectively. We have Re A d = 0.823+°^ i m A rf = -0.199+°^, Re A s = 0.965+g^|, and Im A s = 0.00 ± 0.10. The p-value 
for the 2D SM hypothesis A d = 1 (A, = 1) is 3.0 cr (0.0 cr)'. 



results have changed the constraints on |A,| and our fit is 
significantly worse than in 2010 lfl9l . The value of 0, 
prevents large contributions to Asl from B s so that the 
two measurements do not show a good compatibility. 
We obtain pull values for Asl an d <P^ ~ 2/3. s of 3.0 cr and 
2.7 cr respectively (it was 1.2 cr and 0.5 cr in Ref. |[T9"ll ). 

4. New Physics in A.F = 2 operators: T 12 

In view of these results, new CP -violating contribu- 
tion to rj 2 can be advocated to explain the D0 mea- 
surement of Asl OT1 . It amounts to postulating new B s 
decay channels with large branching fraction. This in- 
duces a SM deviation in AY S as T s n is dominated by the 
CKM-favoured tree-level decay b — > ccs. Any new de- 
cay mode increases the total B s width, whereas the mea- 
sured r s /T d = 0.998 + 0.014±0.012 agrees with the SM 
expectation < r s /T d —l < 4-10~ 4 . The new interaction 
opens b — > s decay modes affecting precisely measured 
inclusive B d and B + quantities [19|. New decays medi- 
ated by a particle with mass M > Mw will add a term 
of 0{M A w /M A ) to ryr™' s , whereas A s normally gets 
a larger contribution of 0(M^,/M 2 ). In models enhanc- 
ing modes with a fermion pair (/, /), one can solve this 
problem with chirality suppression. The extra contribu- 
tion to MJ 2 is down by a factor of mj/M 2 , while that to 
is affected by the milder factor of rrA/mh However 



quantities like Td, s will not be chirality suppressed. It 
seems thus hard to add large NP effects to H, without 
spoiling the agreement of well-measured with the SM. 

Phenomenologically it is much easier to postulate 
NP in T% as it stems from Cabibbo-suppressed decay 
modes like b — > ccsd. Chirality suppression can be in- 
voked to avoid problems with Mf 2 , and there is no prob- 
lem with inclusive decay observables like the semilep- 
tonic branching fraction or the unmeasured AT^. In 
Ref. 1201 , we have studied a new scenario (Scenario IV) 
including the possibility of NP in I^ 2 . It enables NP 
in the |AF| = 1 transitions contributing to Tj 2 , but not 
in other |AF| = 1 quantities entering our fits, such as 
BR(B — > tv). No new CP phase in b — » ccs, which 
would change <p^ , is considered. We introduce new pa- 
rameters 

5 q = ryM^/Re (itf/AO q = d, S (13) 

Re S q , Im S q amount to {AY q l AM q )l{AY™ I AMf 1 ) and 
-flg L /(Ar^ M /AM^ M ), respectively. Re cv is experimen- 
tally only weakly constrained and the correlation be- 
tween Re 5 S , Im 5d and Im 5 S is shown in Fig. [3] The p- 
value of the 8D SM hypothesis Ad = A s = 1, 5 d , s = 
is 2.6 cr. Too large values for \6 S - <Sf M | would be in 
conflict with other observables as explained above. One 
can also study this scenario without NP in the B s sector 
(A, = 1 and 6 S = 5 s ^m)- It could accommodate the main 
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Figure 3: Constraints on (Im Sj, Im <5,) and (Re S s , Im 6 S ) in Scenario IV. The 68% CL intervals are Im = 0.92*1H, Im S s = 1.2*} g. The best-fit 
values of the SM predictions are S S , M = 1 + 0.097 i and Sf* = 1 - 0.0057 i, and the p-value for the 2D SM hypothesis Im 5 d , s 



Im^ M is3.2<T. 

d,s 



anomalies by improving the fit by 3.3<x, but with large 
contributions to Y d n : Im 5d - 1.60*Hg. 

Considering the difficulties to accommodate all in- 
puts (in particular B s — > J/if/<[> and the dimuon charge 
asymmetry, independently of the discrepancy between 
BR(B — > tv) and sin(2/3) H]|6][8)), new measurements 
of flj L , a s Sh and/or a s Sh - a^ h are necessary to determine 
whether scenarios with NP in T d l2 and/or rf, are a vi- 
able explanation for the discrepancies in neutral-meson 
mixing observables with respect to the SM, or if other 
avenues should be investigated (e.g., Ref. Il24l0 . 
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